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We present a numerical study using a cluster algorithm for the 1-d 5 = 1/2 quantum Heisenberg models. 
The dynamical critical exponent for anti-ferromagnetic chains is z = 0.0(1) such that critical slowing down is 
eliminated. 



1. INTRODUCTION 

Since 1976 it is well-known that several map- 
pings of the partition function of a d-dimensional 
S' = 1/2 quantum spin model onto a (d + 1)- 
dimensional classical spin model exist §§]. The 
maps lead to a Feyman path integral formula 
of quantum statistical mechanics. The equiva- 
lence opens the possibility to use powerful com- 
putational techniques to obtain information on 
quantum spin systems. Monte Carlo (MC) sim- 
ulation, a method to get precise numerical data, 
has played an important role in the studies [^, 
and has especially been applied to the anti- 
ferromagnetic models |^J^ since Haldane's con- 
jecture |p[ and the discovery of high- Tc supercon- 
ductors [0 . There are, however, some problems 
with the standard numerical methods, which use 
importance-sampling techniques for the {d + 1)- 
dimensional induced classical spin systems which 
have multi-spin interactions. Their effect is char- 
acterized by large autocorrelation times (in units 
of MC sweeps needed to create a new statistically 
independent configuration) in the low tempera- 
ture region and in the euclidean time continuum 
limit. In this paper a cluster algorithm is devel- 
oped for simulations of quantum spin systems]^, 
especially of 1-d quantum Heisenberg models, af- 
ter mapping them to 2-d induced classical spin 
models and using a blockspin scheme. 



Figure 1. Heisenberg magnet on a chain with an 
even number of spins 2N. 



2. METHOD 

Consider for definiteness a 1-d S* = 1/2 Heisen- 
berg model defined on a chain with an even num- 
ber of spins L = 2N (fig.l) with periodic bound- 
ary conditions. The Hamilton operator is given 
by 



H 



J 



i+i- 



(1) 



Here ai is a Pauli spin operator at the point i 
on the chain. J > corresponds to an anti- 
ferromagnet (AF), while J < corresponds to 
a ferromagnet (F). Using the checkerboard de- 
composition H — Hi + H2 (fig-1) and the Trotter 
formula the partion function is expressed as 
a path integral given by 

Z = Trexp(-/?i/)=[] ^ exp(-5), (2) 

x,t s{x,t) = ±l 
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where (3 is the inverse temperature. Eq.(2) is 
equivalent to a 2-d classical spin model on a rect- 
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angular lattice L x 2M with Ising-like variables 
s{x,t) — ±1. The lattice spacing in the ad- 
ditional so-called euclidean time or Trotter di- 
rection is e = 1/M (fig. 2). The classical spins 
interact with each other via 4-spin couplings 
exp(-S'[si, S2, S3, S4]) = (si, S2I exp(-e/3J(7i • 
(Ti+i/4)ls3, S4), which are given by elements of the 
transfer matrix. For the AF one has 

S[l, 1, 1, 1] = -1, -1, -1] = e/3J/4, 

S-il, -1,1,-1] = S'[-l, 1,-1,1] = 

e/3J/4-log[(exp(e/?J) + l)/2], 
S-il, -1,-1,1] = S'hl, 1,1,-1] = 

e/3J/4-log[(exp(e/3J)-l)/2]. (3) 

All other action values are infinite and the cor- 
responding elements of the transfer matrix are 
zero. The matrix elements can be interpreted as 
the Boltzmann weights of the spin configurations. 
Up to here the quantum spins Ui with 2-spin cou- 
plings living on a 1-d chain have been mapped to 
classical spins s{x,t) with 4-spin couplings living 
in the 2-d plane. However, most classical spin 
configurations are forbidden, because for them 
some elements of the transfer matrix are zero 
(their Boltzmann factor vanishes). Therefore, it 
is natural to attempt to choose a collective non- 
local update technique. For models with 2-spin 
couplings this can be done using the Swendsen- 
Wang 1^ or Wolff jl^] cluster algorithms which 
flip whole cluster of spins simultaneously. These 
algorithms, however, can not be applied directly 
to a model with 4-spin couplings. 

To make an application of the cluster tech- 
nique possible we further map the classical spin 
model with 4-spin couplings S'[si, S2, S3, S4] to 
a blockspin model. A blockspin B{n, m) ~ 
{s{x,t),s{x + l,t),s{x,t + l),s(x + l,t + 1)} is 
a collection of four spins, where x = 2n — l,t = 
2m. The blockspins interact via 2-blockspin in- 
teractions S[B,B'] — S'[si, S2, S3, S4] as shown in 
flg.2a. Each spin belongs exactly to one block- 
spin and the blockspins live on a lattice with a 
doubled lattice spacing. The spins can also be ar- 
ranged to blockspins in another way, B(n,m) = 
{s{x , t) , s{x + 1, t) , s{x , t + 1) , s{x + 1, t + 1)} with 
X — 2n, t — 2m — 1 (see fig. 2b). An updating algo- 
rithm for the blockspins must alternate between 



Figure 2. 2-d induced classical spin system with 
4-spin couplings depicted by shaded plaquettes. 
(a) The scheme {B}: each blockspin consists of 
4 spins and interacts with its nearest-neighbors 
via the original 4-spin interaction, (b) the other 
scheme {B} which is dual to the scheme {-B}. 



the schemes {B} and {_B} to ensure ergodicity. 
We have also chosen two kinds of nonlocal block- 
spins to change the magnetization and the so- 
called winding number by odd integers. The clus- 
ter algorithm makes use of a flip symmetry of the 
blockspin model. A blockspin B — {si, S2, S3, S4} 
is flipped to —B = {— si, — S2, — S3, — S4} simply 
by flipping all spins in B. The blockspin action 
B'] is invariant under a flip of all blockspins 
because the original action is symmetric against 
flipping all spins. 

3. MONTE CARLO RESULTS 

We now describe how to generate transitions 
between blockspin conflgurations. a) The algo- 
rithm puts bonds between all nearest-neighbor 
pairs of blockspins with a probability p = 
1 - min{l, exp(-5[-B, B'])/ exp(-5[B, B'])}. b) 
Two blockspins belong to the same cluster if they 
are connected by a bond, c) All blockspins in one 
cluster must be flipped simultaneously, d) Us- 
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Figure 3. A log-log plot of the autocorrelation 
time Te versus 1/e at (3 J — 1 both for (a) AF and 
(&) F Heisenberg magnets with 2N = 32. 



ing the Swendsen-Wang multi-cluster method [£| 
each cluster is flipped independently with proba- 
bility 1/2. For the Wolff single-cluster algorithm 
|l0| one blockspin is randomly chosen and the 
cluster to which it belongs is flipped. The pro- 
cedure a)-d) of changing a configuration satisfies 
the detailed balance condition ||l^ . For ergodicity 
we alternate between the {B} and {B} schemes. 
We tested the blockspin cluster algorithm in de- 
tail for one-dimensional spin chains. To demon- 
strate the efficiency of the algorithm we compare 
it to a blockspin Metropolis update. For both 



algorithms we measure the autocorrelation func- 
tions Ca{t), where O denotes the observables, i.e. 
the internal energy e, and the uniform and stag- 
gered susceptibilities x ^^'^ Xs- Then we obtain 
the integrated autocorrelation times To by 

oo oc 

eM~l/ro)=J2^o{t)/Y.Co{t). (4) 

i=l t=Q 

In all the simulations we have performed a ran- 
dom start followed by 5000 sweeps for thermaliza- 
tion and 50000 sweeps for measurements using the 
single-cluster updating. From the data at /? J = 1 
shown in fig. 3 we obtain the dynamical critical ex- 
ponent Ze both for the anti-ferromagnet and for 
the ferromagnet. z^ is defined in the continuum 
limit e ^ by Te oc We find 

Ze = 0.0(1) for the cluster algorithm, 

Ze = 0.8(1) for the Metropolis algorithm. (5) 

The autocorrelation times t-^ and t^^ do not 
diverge in the continuum limit even for the 
Metropolis algorithm. Still, the corresponding 
autocorrelation times are at least an order of mag- 
nitude larger than the ones of the cluster algo- 
rithm. 

As shown in fig. 4, the Metropolis algorithm has 
severe problems with slowing down at low tem- 
peratures. At f3J = 8 for example t^^ = 12(1) 
for J > 0, and = 3300(200) for J < 0. The 
cluster algorithm, on the other hand, has auto- 
correlation times of at most a few sweeps and it 
can update the configurations efficiently at lower 
temperatures like e.g. f3J = 16. In fig. 4 the 
and are fitted by 

oc exp(cx/3|J|), T^^ ocexp{c^J\J\). (6) 

The data at /3 J = 16 show some slowing down for 
the ferromagnet, but the autocorrelation times 
are moderate. For anti-ferromagnet there is no 
indication of slowing down such that 

2^=^^=2^^=0.0(1). (7) 

4. DISCUSSION 

First, we like to mention a few words about 
frustration. As it is well known the Swendsen- 
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Figure 4. (a) The plot of logr,^^ versus (3 J for 
AF couplings and {h) logr^ for F couplings with 
2N = 128 and e = /3/M 0.125. 



Wang and Wolff cluster algorithms eliminate crit- 
ical slowing down for the Ising and Potts mod- 
els because they are not frustrated. For strongly 
frustrated models, e.g. for spin glasses, cluster 
algorithms do not work efficiently. Therefore the 
question arises if our blockspin models are frus- 
trated or not. We have verified that the blockspin 
model for 1-d AF couplings has indeed no frus- 
trated allowed configurations. For the 1-d ferro- 
magnet, on the other hand, some allowed config- 
urations have a weak frustration. However, in the 
continuum limit e — > the frustration disappears. 



We believe that the weak frustration is responsi- 
ble for the mild slowing down in the ferromagnetic 
case at low temperature. 

The blockspin cluster algorithm can also be 
applied to other models, e.g. to systems with 
anisotropic couplings or in an external magnetic 
field, to models of higher spins (s = 1, 3/2, ...), to 
six- and eight-vertex models and to 1-d fermion 
systems. Work on the 2-d anti-ferromagnet is in 
progress [ pl| . 
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